Local phase damping of single qubits sets an upper bound on the phase damping rate 

of entangled states 
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I derive an inequality in which the phase damping rates of single qubits set an upper bound for the 
phase damping rate of entangled states of many qubits. The derivation is based on two assumptions, 
first, that the phase damping can be described by a dissipator in Lindblad form and, second, that 
the phase damping preserves the population of qubit states in a given basis. 
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I. INTRODUCTION 



Quantum information processing [l| offers a perspec- 
tive for a tremendous reduction of the computation time 
in solving certain problems, such as factorization of large 
numbers [2], simulations of quantum systems Q, and 
database searches However, interactions of the quan- 
tum system with its environment induce decoherence [B- 
[^. This is a major limiting factor on the way toward 
large-scale experimental implementations. Naively, one 
might expect that the decoherence rate of an entangled 
many-qubit state should equal the sum of the decoher- 
ence rates of the individual qubits. While this is true if 
the decoherence processes are local, i.e. if they act inde- 
pendently on individual qubits, the situation can change 
if the decoherence processes act on some or all qubits in 
a correlated manner. Such correlated decoherence can 
give rise to decoherence-free subspaces which have been 
studied theoretically [ol-fTsI and experimentally fl3 - [l7j . 
This shows that certain entangled states can decohere 
considerably slower than their single-qubit constituents. 
In other words, a single-qubit decoherence rate is not a 
lower bound for the decoherence rate of entangled states. 



Here I show that the single-qubit decoherence rates 
set a rigorous upper bound for the decoherence rate of 
entangled states of n qubits. The existence of such an 
upper bound is somewhat surprising because entangle- 
ment is difficult to prepare and maintain experimentally. 
The upper bound derived here is experimentally relevant 
because in many experiments where one aims at gener- 
ating entangled states one has the capability to measure 
single-qubit decoherence rates already during the build- 
up phase of the apparatus, often much earlier than the 
time where one manages to generate and detect entangled 
states. In addition, one often uses experimental tech- 
niques such as a magnetic hold field to suppress spin-flip 
transitions of the qubits. But these techniques cannot 
protect the qubit against loss of phase coherence. The 
experimenter thus selects a preferred basis in which the 
populations are preserved, whereas phase coherence may 
decay. 



II. MODEL 

The derivation of the upper bound rests on two as- 
sumptions. The first assumption is that there is an or- 
thonormal basis B in which the population of the basis 
states is time independent whereas the relative phases 
between the basis states may rotate and decay. This 
condition is easy to verify experimentally. I will show 
below that this requires the Hamiltonian H to be diago- 
nal in the basis B. This implies that the final results will 
be relevant for quantum memories, but they are not nec- 
essarily applicable during quantum gate operations that 
act on the qubits of interest. The second assumption is 
that the decoherence can be described by a dissipator in 
Lindblad form. This is a reasonable assumption for many 
experiments, as discussed now. 

Consider a quantum system coupled to an environ- 
ment. If one assumes, first, that the system is initially 
not entangled with the environment and, second, that 
a Born-Markov approximation is appropriate (because 
the coupling between system and environment is weak 
enough and because the environment is much larger than 
the system), then the environmentally- induced decoher- 
ence can be described by 
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[H,p]+Vp 



(la) 



with a dissipator V in Lindblad form |18M2' 



_,(rn) , 



.(^(m))t^(m)^_^(^(m))t^(m)^ . (lb) 



Here, p is the density matrix, the 7^'") > are deco- 
herence rates, and the dimensionless operators yl^™) are 
called Lindblad operators. In the following, it is always 
assumed that the time evolution is described by Eq. ([l}. 

Definition: I call a time evolution generated by Eq. 
([1]) population preserving with respect to an orthonormal 
basis B if all diagonal elements of p in the basis B are time 
independent for all initial states, i.e. {d/dt)pn = 0. 
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III. RESULTS 

I now formulate three theorems, which are proven in 
appendix [Xj 

Theorem 1 : If aU the yl'™) are diagonal in the same 
basis B with eigenvalues A , then 

Vp,j ^ {iA,j -rij)p,j (2) 



with decay coefficients Tij and angular frequencies 
given by 



,(™) 



^(m) ^(m) 



(3) 



(4) 



If additionally, H is diagonal in the same basis with eigen- 
values Ej, then 



dt 



with angular frequencies 



E-i — Ei 



(5) 



(6) 



Note that > and Tu = An ~ uju ~ 0. Hence, only 
off-diagonal elements of p change over time. These ele- 
ments do not mix. Instead, each off-diagonal element ex- 
periences two effects. The first effect, described by Tij, is 
an exponential decay of \pij\, which is called phase damp- 
ing The second effect, described by w^j, is a phase 
rotation [2l[. This paper focuses on phase damping. 

Theorem 2: A time evolution is population preserving 
with respect to B, if and only if H and all the A*^™^ are 
diagonal in B. 

Theorem 5: If a time evolution is population preserv- 
ing with respect to B, then the Tij from theorem 1 obey 
the following inequality for n e {1,2,3,...} and for all 
basis states \po) e B, \pi) e S, . . . , |p„) e B 



■Pk- 



(7) 



k=l 



This inequality becomes an equality, if and only if for 



each m, the expression Ap™jj — Ap"''' is independent of fe, 
where the Aj-™' are the eigenvalues of the A^'^\ 

This theorem is applicable to the density matrix p of 
an arbitrary state of n qubits. Let the states | f) and | \) 
denote a basis of the Hilbert space of each single qubit 
and assume that the time evolution is population pre- 
serving with respect to the basis B which consists of all 
the tensor products of the single-qubit states | f) and | ],) . 
This situation corresponds, e.g., to the scenario with a 
magnetic hold field discussed in the introduction. Ac- 
cording to theorem 1 , each density matrix element pij in 
the preferred basis B obeys Eq. ([S]). Specifically, there 



(m) 



will be phase damping, described by Tij. For any choice 
of the basis states |i) and \ j) one can obviously construct 
a sequence of basis states jjio) , \pi), ■ ■ ■ , \Pn) in a way that 
bo) = \Pn) = and that (for all fc G {1, 2, . . . , n}) 
the states \pk-i) and \pk) differ by only one local spin 
flip between states 1 1) and 1 1) . Application of Eq. ([7]) 
yields an upper bound for Tij and each Tp^^ -^^p^, on the 
righthand side obviously describes a local phase damping 
rate. 

This is the central result of the present paper. Interest- 
ingly, the upper bound in the inequality (O is a factor of 
n higher than the decoherence rate that one would obtain 
for uncorrelated phase damping of individual qubits. 

The experimental relevance of this upper bound arises 
from the fact that measurements of the single-qubit phase 
damping rates are fairly easy to perform because no en- 
tangled state needs to be prepared and detected. Accord- 
ing to the inequality such measurements already set 
a worst-case upper bound for the phase damping rate of 
any entangled state. 



IV. EXAMPLES 

To illustrate this concept, consider an example of a 
Greenberger-Horne-Zeilinger(GIIZ)-type entangled state 
[2^ of n qubits 



|V'GHz) = ^(|t)®" + l;)^ 



(8) 



Here, | -1^) ®" = | f) ® • • • (g) | f) abbreviates a tensor product 
of n times the same single-qubit state and the inequality 
([7]) can be applied to the basis states 



\Pk) = |t)®fe® 



(9) 



with k G {0, 1, . . . , n}. This yields an upper bound for 
the phase damping rate of the GHZ state because Fp^.p^ 
describes the decay of {po\p\Pn) = (i ■ ■ • i |p| t ' ' ' t)- 
The other coefficients Fp^ ^ .p^ appearing in the inequal- 
ity ([7]) describe decay of {pk-i\p\pk}- The states \pk-i) 
and \pk) differ only by a spin flip of the k-th qubit. Hence, 
this coefficient can be determined experimentally from a 
measurement of the single-qubit phase damping rate of 
the fc-th qubit. 

Let Fghz = Tpo.pn denote the phase damping rate 
coefficient of the GHZ state and let F^ = Fp^ ^^.p^ de- 
note the phase damping rate coefficient of the fc-th qubit, 
where the specific orientation of the other qubits was 
dropped from the notation for brevity. Then Eq. ([7]) 
yields 



Fghz < n(Ti + F2 



F„). 



(10) 



To illustrate this inequality further, consider some ex- 
amples for n = 2 qubits. Here, one typically denotes the 
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Bell states as 



^/2 
1 

71 



(in>±iu». 



(11) 

(12) 



For 2 qubits, the GHZ-type state of Eq. ([5]) is obviously 
the Bell state 

Example 1 : Consider local phase damping generated 
by the Lindblad operators 



^^'^ = It)(t|l®l2, 
= li®|t>(t|2 



(13) 
(14) 



where 1 denotes the identity matrix. This yields Fi = 
7(1), F2 =7(2), and 



(15) 



Local phase damping acts identically on all Bell states. 
The phase damping rate of each Bell state is simply the 
sum of the local phase damping rates. 

Example 2 : Consider phase damping generated by the 
Lindblad operator 



^^'^ = in)(ni-i;;)(;;i 



(16) 



The states |tt) and 114-) experience phase damping rel- 
ative to the two-dimensional rest of the Hilbert space. 
Here, Fi = F2 = t^^^ and 



0, 



F.± =2(Fi-HF2). 



(17) 



The states IV'^) span a decoherence-free subspace, 
whereas the states experience phase damping at a 
rate that reaches the upper bound ([T0|) for n = 2. 

This example is somewhat related to the experiment 
with two trapped ions in Ref. 15[. In that experiment, 
the phase damping is dominated by fluctuating ambient 
magnetic fields with frequencies primarily at 60 Hz and 
its harmonics. These fields cause a fiuctuating Zeeman 
energy for the ions. The ions are separated by only a 
few micrometers so that the ambient magnetic fields are 
roughly uniform across the trapping region. As a result, 
the states | fi) and | ]/[) experience no net Zeeman effect 
and the states \tp'^) experience no phase damping, at least 
to lowest order. The states |tt) and |||), however, expe- 
rience plus or minus twice the single-qubit Zeeman shift, 
leading to a phase damping that is twice as fast as for a 
single qubit. 

Note that technical fluctuations of a macroscopic mag- 
netic field need not necessarily allow for a description in 
terms of Eq. ([T]). However, if a similar experiment were 
performed with quantum dots in a solid, then the mag- 
netic dipole-dipole interactions with the large number of 
surrounding nuclear spins with thermal occupation may 
create an environment that can be described by Eq. ([l}. 

Example 3: For comparison, consider phase damping 
generated by the Lindblad operators 



AW = |tt)(ttl, = !;;)(;; I 



(18) 



with 7 



(1) 



.(2) 



Here, A^^^ creates phase damping for 



the state | ft) relative to the three-dimensional rest of the 
Hilbert space and j4(2) creates an analogous effect for the 
state lU)- This yields Fi = F2 = 7^^^ and 



0, 



F1+F2. 



(19) 



Again, the states span a decoherence-free subspace. 
But now, the phase damping rate of the states \<j)^) no 
longer reaches the upper bound (fTO|) . 

Such a phase damping might be obtained hypotheti- 
cally, if there are two different fields, one coupling only 
to state I tt) and another coupling only to state 1 1|) . If 
these fields are uncorrelated, they will each generate an 
individual Lindblad operator, as in Eq. (jl8p . This situa- 
tion differs from the phase damping caused by a common 
field, which is expressed by one Lindblad operator as in 

Eq. (nr 
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Appendix A: Proofs 



Proof of theorem 1 : A 



(m) 



yield Pp., = '-^S,,pu6i,{2Xt\x^"'^ 

v(m)|2 



and Eq. ([T]) 

(™) 12 



I )i where Sij denotes the Kronecker symbol. Com- 
bination with |A,("' - Aj™^^ = |a|™)|2 + lA^-^'p - 

2Re(A|"^(Aj™')*) yields Eq. @. Combination with 
{d/dt — 'D)pij = {EiPij — pijEj)/ih completes the proof. 

Proof of theorem 2: li H and all the A^™) are diagonal 
in B, then theorem 1 yields F^^ = uja = and {d/dt)pii = 

for all i. 

To prove the other direction, assume that the time 
evolution is population preserving with respect to B. I 
will first show that all the A^"^'^ are diagonal in B. To 
this end, it suffices to consider {d/dt)pjj = for the 
special case where p = \j){j\ at i = where \j) is any 
pure basis state that belongs to B. Eq. (P) yields — 

{d/dt)p,, = - 7^") Efe^, l47^ I'- This holds for all 
j. This yields ^i™'' = for all j, k, m with j 7^ k. Hence, 
all the A*^™^ are diagonal in B. 

It remains to be shown that H is diagonal in B. To 
this end, select an arbitrary pair of indices i,j with 

1 j- It needs to be shown that Hij = 0. Here, it 
suffices to consider only initial density matrices of the 
form p = with = + e''^\j))/V2 with a 
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real parameter i^. As seen above, all the ^('"^ are diag- 
onal in B so that theorem 1 yields Vpa — 0. The pop- 
ulation preserving character of the time evolution then 
yields = ifi{d/dt)pii = Y.ke{i^j}iHikPki - PikHki) = 
HijPji — c.c. = ilm{Hije'^'-P). This holds for all real (p, 
especially for p = and p — tt/2. Hence Hij = 0. 

Proof of theorem 3: Assume that the time evolution 
is population preserving. According to theorem 2, all 
the are diagonal in B. Hence, theorem 1 is appli- 

cable. The Cauchy-Schwarz inequality \ Y^k=i ^kUll'^ < 



(ELil^fcP)(ELilyfcP) with xk 



Vk 



IreadslA^r^-ACf = IELi(A; 



(™)|2 



Am) 



, (m) 
>k-l 
. (m) 



(m) 



and 



Aj,™Y- Multiply this by 7('"V2- Then 
sum over m to obtain the inequality ([T]) ■ It is well known 
that the Cauchy-Schwarz inequality becomes an equality, 
if and only if the vectors (xi , . . . , a;„) and , . . . , 7/„) are 
linearly dependent. Here, this is equivalent to xi — X2 — 
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